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Abstract 

We consider stochastic diffusion processes absorbed at the boundary of a domain. 
It is shown that there exist initial distributions which ensure a given decreasing of 
density of the absorbed process. 
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1 Introduction 

The paper investigates distributions of diffusion processes with adsorption on the boundary 
of a domain. We consider the problem of controlled absorption, and we obtain a new 
property of absorbed process: we establish an existence of initial distributions which 
ensure a given decreasing of density of the absorbed process. More precisely, we show that, 
for any piecewise continuous function 7(2;) > such that 7Q 7^ 0, there exists an initial 
probability density function 0) and a number a > such that p(x, 0) = p(x, T)+a r y(x), 
where p(x, t) is the probability density function of a process y(t) absorbed at the boundary 
dD. Thus, it can be concluded that decreasing of density of the absorbed process can be 
programmed. 

It can be mentioned that some problems for controlled absorption were solved by author 
(1994, 1995) for another type of constraints, when v(x, T) = fiv(x, 0) with /i € (0, 1), where 
v(x,t) is a solution of a backward Kolmogorov equation. 
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2 Definitions 

Consider a n-dimensional Wiener process w(t) with independent components on a complete 
probability space (fi, J 7 , P), such that io(0) = 0. 

Let D G R n be a bounded domain with C 2 - smooth boundary dD, and let T > be 
a fixed number. Consider the following Ito's stochastic differential equation: 



Here yo is a random vector which is independent of w(-) and such that y € D a.s. 

We assume that the function f(x,t) : R n x R — > R n is measurable and bounded, the 
function 0(x,t) : R n x R — > R nxn is continuous, and that there exist bounded derivatives 
d k (5{x, t)/dx k , where i = l,...,n, k = 1,2. We assume also that (3(x,t)P(x,t) T > 5I n , 
where 5 > is a constant, and I n is the unit matrix. 

Under these assumptions, there exists the unique weak solution y(t) of (|l|) (see e.g. 
Gihman and Skorohod (1979)). We consider the process y(t) which is absorbed at dD, 
i.e. until first exit from D. 

Spaces and classes of functions 

For a Banach space X, we denote the norm by || • \\x- 

o 

Let H° = L 2 {D) and H 1 =W\ (D) be the standard Sobolev Hilbert spaces. Let H' 1 
be the dual space to H , with the norm || • ||#-i such that for u G H° is supremum 

of (u,v) H o over all v G H° such that ||v||#i < 1; we have H 1 C H° C H , assuming the 
standard embedding. 

We shall denote the Lebesgue measure and the <7-algebra of Lebesgue sets in R m by 
l m and B m , respectively. 
Introduce the spaces 



dy(t) = f(y(t), t)dt + (3{y(t), t)dw{t), t E [0, T] 
2/(0) = yo- 



(1) 



C k (s,T)^c([s,T];H k ), W k (s,T)^L 2 ([s,T],B 1 ,li;H k ), k = 0,±l 



and the space 



V k {s,T) = W k (s,T)nC k " 1 (s,T), k>0 



with the norm ||tt||v*( s ,T) = ll n llw fe (s,T) + ||ti||cfc-i ( s ,t) - 
The following definition will be useful. 
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Definition 1 A function j(-) : D — > R is said to be piecewise continuous if there exists a 
integer N > and a set of open domains {Di}f =l such that the following holds: 

• D = uf =1 D u Di n Dj = for i ± j; 

• for any i € {1,...,^}, the function 7^ is continuous and can be continued to a 
continuous function 7^ : Di U dDi — > R; 

• for any x E ufL^dDi, there exists j € {1, ...,iV} snc/i i/jai x € <9.Dj and 7 jj(x) = 7(2;). 



3 The result 

Theorem 1 Ze£ 7(-) : R n — > R 6e a piecewise continuous function such that j(x) > 
(Vx) and 7(-) 7^ 0. T/ien i/tere exists a random vector y$ which is independent ofw(-) and 
such that the following holds: 

(i) the process y(t), defined by ([7p and absorbed at dD, has the probability density 
function p{x,t) £ V 1 (0,T); and 

(ii) iaere exists a constant a > swca i/iai 0) = p(x, T) + 07(2;). 

Proof. Let yo have a probability density function p(x) € L2(D). In that case, the 
probability density function p(x, t) of the process y(t) absorbed at dD satisfies the the 
forward Kolmogorov equation 



(2) 



% = Ap, t>0, 
p(x,t)\ xEdD = 0, p(-,0) = p(-), 

where 

n q2 n d 

A P=Y, dx . dx . Mx,t)p(x))-^-Q-(fi(x,t)p(x)) 

i,j=l i 3 i=l ' 

and a(x,t) = {aij(x,t)} = \(3(x, t)(3{x, t) T . 

For s G [0, T] and ^ E consider the following auxiliary boundary value problem: 

S = ^> (3) 

«(s,i)UeaD = 0, v(-,s) =£(•). 
From the classic theory of parabolic equations, it follows that 

\\v(;T)\\ m KCtM-^^Hl, (4) 
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where C\ > is a constant which does not depend on £ and s (see e.g. Ladyzenskaya ei 
al (1968)). 

Introduce operators £ s : H° — > V 1 (s,T), such that £ s £ = w, where u is the solution in 
V 1 (s,T) of the problem (Q). These linear operators are continuous (see e.g. Ladyzenskaya 
et al (1968)). Introduce an operator Q : H° — ► if , such that Q£ = v(-,T), where v = £o£- 
Clearly, this operator is linear and continuous. 

Lemma 1 (i) The operator Q : H° — > H° is compact; 

(ii) If the equation Q£ = £ has the only solution £ = in H , then the operator 
(I — Q) _1 : H° — > H° is continuous. 

Proof of Lemma [J. Let £ € H° and v = £o£> i- e - u is the solution of the problem (|3|). 
We have that u| tg [ Sj T] = C s v(-,s) for all s € [0, T], hence 



K-,T)|| H i < Ciinf te[0iT ] \\v(-,t)\\ H i 

% {lo H;t)\\ 2 mdt) 1/2 < VL\\v\\ wl < ^LU\\h° 



< 



for constants Cj > which do not depend on £. Hence the operator Q : ii° — ► H 1 is 
continuous. The embedding H 1 to if is a compact operator (see e.g. Yosida (1965), 
Ch.10.3). Then (i) follows. Further, (ii) follows from Fredholm Theorem. This completes 
the proof of Lemma [l]. □ 

Lemma 2 For any 7 € H°, there exists the unique solution u € V 1 of the following 
problem: 

' 8 = An, t > 0, (s) 
«(a,*)|xe9D = °> «(x,0) - -u(x,T) = 7(2;). 

Proof of Lemma\^. First, we shall show that if 'y(-) = then the unique solution of (|5|) 
in V 1 is u(-) = 0. 

Let u G V 1 solve © with 7(-) = 0. Clearly, u = C Q u(-,0). Denote C + (^) = 
max(0, u(x, 0)) and C~(a;) = max(0, — u(x, 0)). Denote ii~ = £oC~ an d ' u+ = £oC + - 
We have that u + > and u _ > a.e., u = u + — u _ and u(x, 0) = ( + (x) — £~(x). 

If «(-,0) 7^ then either £ + (-,0) 7^ or £~ (-,0) 7^ 0. It follows from absorption at dD 
that if C + (-,0) / then 

/ u + (x,T)dx < / ( + (x)dx. 
Jd Jd 

(It suffices to note that the process q = v^u^ is a probability density function of a process 
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with absorption on 5D, where v = J D ( + (x)dx). Similarly, if ( (-,0) ^ then 

/ u~(x,T)dx< / Q~(x)dx. 

Hence 

J D \u(x,T)\dx < j D u + (x,T)dx + f D u~ (x,T)dx 

< f D C + {x)dx + f D C~(x)dx = j D \u(x,0)\dx, 

i.e. f D \u(x, T)\dx < J D \u(x, Q)\dx, and the condition u(x, 0) = u(x, T) fails to be satisfied 
for u(-) ^ 0. Thus, the unique solution of (g) for j(-) = is u(-) = 0. 

By Lemma [l], it follows that the operator (/ — Q) _1 : H° — > i7° is continuous. Let 
7(0 G H° be arbitrary. Then there exists £ = (J — Q) _1 7 € -fT , and this £ is unique. Let 
u = £o(. By the definitions of Co and Q, it follows that u(-,T) = Qu(-,0). We have that 
u(; 0) - Qm(., 0) = t(-), i-e. «(-, 0) - «(-, T) = 7 (0- Thus, u = £ ( = £ Q (I - Q) -1 7 is the 
unique solution of (JsJ) for any 7(0 € H° = L2(D). This completes the proof of Lemma ||. 
□ 

Let us continue the proof of the theorem. Let j(x) > be such as in the assumptions 
of the theorem, and let u = C${I — Q) _1 7 solve the problem (|5|). 

It is easy to see that if u(-, 0) = then u(-,T) = and 7(0 = 0. By the assumpltions, 
7(0 + 0, hence «(-,0) + and «(0 + 0. 

We remind that u = CoC, where C = u(-,0) E H°. By Theorem 9.1 from Ch.IV of 
Ladyzenskaya et al (1968) applied for smooth functions which approximate u(-,0) in H°, 
and by Theorem 8.1 from Ch.III of this book, it follows that there exists a representative 
u(-,T) of the corresponding element of H° which is continuous in x G D, where D = 
D U dD. (Note that, by the definition, an element of H° = L2(D) is a class of £ n - 
equivalent functions). We have that u(-,0) = u(-,T) +7(0, hence there exists a piecewise 
continuous representative u'(-,0) of u(-,0) € H°. Note that such u'(-,0) is not unique, 
because the boundary values at dDi in the Definition [l] can be choosen differently. 

Let us show that u(x, 0) > for a.e. x. Suppose that 

3x G D : u'(x,0)<0. (6) 

If (||) holds, then there exists a piecewise continuous representative u(-, 0) such that there 
exists x G D such that 

u(x, 0) < 0, u(x, 0) < «(x, 0) for a.e. x £ D. 
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We have that 

u(x,T) = G(x,y,T,0)u(y,0)dy, 
Jd 

where G(x, y, T,0) > is the corresponding Green's function for the problem (jl]). Clearly, 

/ G(x, y, T, 0)dy <E (0, 1). 
Jd 

Then 

u(x,T) = (f D G(x,y,T,0)dy^u(x,0) + J D G(x,y,T,Q)(u(y,0) -u(x,0))dy 
> ru(x,0) > u(x,0). 

It follows that if @ holds then u(x,t) does not satisfy (||). Thus, u(x,0) > a.e. 

Now, let a = (J D u(x,0)dx)~ 1 and p(x) = au(x,0). We have that p{x) > and 
J D p(x)dx = 1. Then there exists a random vector yo such that yo is independent of w(-) 
and has the probability density function p. Let p = Cqp, i.e. p solve (|2|). Clearly, p(x,t) 
is the probability density function of the process y(t) satisfying ([[]). By linearity of (||), it 
follows that p = au and p(-, 0) —p(-,T) = a^f(-). This completes the proof of Theorem |l]. □ 
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